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t^J- ■ ABSTRACT 

Cn| \ We have conducted a programme to determine the fundamental parameters of a sub- 

stantial number of eclipsing binaries of spectral types and B in the Small Magellanic 
Cloud (SMC). New spectroscopic data, obtained with the two-degree-field (2dF) multi- 
object spectrograph on the 3.9-m Anglo- Australian Telescope, have been used in con- 
junction with photometry from the Optical Gravitational Lens Experiment (OGLE-II) 
database of SMC eclipsing binaries. Previously we reported results for 10 systems; in 
this second and concluding paper we present spectral types, masses, radii, tempera- 
tures, surface gravities and luminosities for the components of a further 40 binaries. 
The uncertainties are typically ±10% on masses, ±4% on radii and ±0.07 on logL. 
. The full sample of 50 OB-type eclipsing systems is the largest single set of fundamental 

parameters determined for high-mass binaries in any galaxy. We find that 21 of the 
systems studied are in detached configurations, 28 are in semi-detached post-mass- 
transfer states, and one is a contact binary. 

The overall properties of the detached systems are consistent with theoretical 
models for the evolution of single stars with SMC metal abundances (Z ~ 0.004); 
in particular, observed and evolutionary masses are in excellent agreement. Although 
there are no directly applicable published models, the overall properties of the semi- 
detached systems are consistent with them being in the slow phase of mass transfer 
in case A. About 40% of these semi-detached systems show photometric evidence of 
rS | orbital-phase-dependent absorption by a gas stream falling from the inner Lagrangian 

^ . point on the secondary star towards the primary star. This sample demonstrates that 

case-A mass transfer is a common occurrence amongst high-mass binaries with initial 
orbital periods P < 5 days, and that this slow phase has a comparable duration to the 
detached phase preceding it. 

Each system provides a primary distance indicator. We find a mean distance 
modulus to the SMC of 18.91 ± 0.03 ± 0.1 (internal and external uncertainties; D = 
60.6 ± 1.0 kpc). This value represents one of the most precise available determinations 
of the distance to the SMC. 

Key words: binaries:eclipsing - binaries:spectroscopic - binaries:general - bina- 
ries:close - Magellanic Clouds 
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1 INTRODUCTION 

In a previous paper, Harries, Hilditch & Howarth (2003, 
hereinafter H 3 03) reported the first results of an observa- 
tional programme to determine the fundamental parameters 
of stars in eclipsing binaries in the Small Magellanic Cloud 
(SMC). The rationale for this programme was fully de- 
scribed in H 3 03, so it suffices to state here that the twin 
aims of this investigation are to provide tests of evolutionary 
models for high-mass binary stars, including those that un- 
dergo mass-exchanging interactions at early stages in their 



evolution; and to provide a set of primary distance indica- 
tors to determine the distance to, and line-of-sight depth 
of, the SMC, independently of all other distance calibrators. 
Since publication of H 3 03, a Joint Discussion on 'Extragalac- 
tic Binaries' was held at the 2003 IAU General Assembly, 
and the papers presented there by many authors provide a 
broad overview of recent progress in this field (cf. Ribas & 
Gimenez 2004). 

H 3 03 reported combined spectroscopic and photomet- 
ric analyses of 10 relatively bright binaries; we now present 
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follow-up observations and analyses of 40 additional targets. 
With such a large sample we are insensitive to anomalous 
individual results (whether they arise from astrophysical ef- 
fects or from problems with data, such as unresolved blend- 
ing in the photometry), as well as gaining from simple sta- 
tistical leverage that gives our results (and particularly our 
distance determination) very high weight. Furthermore, by 
introducing radial- velocity data, albeit of modest quality, we 
gain enormously over analyses based solely on light-curve so- 
lutions. 

We summarize our observing programme in Section |2j 
and the spectroscopic and light-curve analyses in Section 
SectionQ]describes comparisons between the derived funda- 
mental astrophysical parameters and evolutionary models 
for single stars and interacting binary stars, while Section 
reviews distance determinations. 



2 SPECTROSCOPIC OBSERVATIONS 

Our programme was devised in 2000, aiming to exploit both 
the publication of the Optical Gravitational Lens Exper- 
iment photometric survey (OGLE-II; Udalski et al. 1998) 
and the commissioning of the two-degree-field (2dF) multi- 
object spectrograph on the 3.9-m Anglo- Australian Tele- 
scope (A AT). Extensive simulations for realistic expecta- 
tions of time allocations showed that satisfactory orbital so- 
lutions should be obtainable for systems with B < 16 mag 
and orbital periods P < 5 days - constraints ensuring both 
adequate phase coverage for a reasonable sample, and suf- 
ficient signal-to-noise ratio on the faintest targets (target 
S/N> 25 per wavelength sample per observation). With 
these boundary conditions (and some further limitations in- 
troduced by practical constraints on the positioning of the 
2dF optical fibres) we identified an initial sample of 169 SMC 
eclipsing binaries. 

Our spectroscopy of these targets was obtained in two 
overlapping 2°-diameter fields in the central regions of the 
SMC; details of the field centres, and checks on the astrome- 
try, are described in H 3 03. The data were taken using 2dF, in 
service mode (as is standard for this instrument). Each full 
observation consisted of two 1800-s integrations (on, typ- 
ically, 100 targets simultaneously), plus arc and flat-field 
frames. Integration times of 1800s were used for the target 
exposures, in seeing of typically 1.5-2" (cp. fibre diameters 
of 2"). We used the marginally better-performing of the two 
spectrographs, with a 1200B grating providing a mean re- 
ciprocal dispersion of l.lApx" 1 and a resolution of 2A (the 
best available with this instrument). The wavelength range 
of the data is 3855-4910A, which ensures excellent cover- 
age of the H, He I and He n lines that dominate the blue 
spectra of O and early B stars. The 2dFDR package (Lewis 
et al. 2002) was used with default settings to reduce all the 
data. The resultant spectra have S/N values ranging from 
~ 20 to ~ 50. 

From allocations of 12 nights over the 2001-2003 SMC 
observing seasons, worthwhile observations were secured on 
only 5j nights. In consequence, the orbital phase coverage 
is generally poorer than planned, and for any given target 
depends on the orbital ephemeris and location (in the field 
overlap region or otherwise). The final database has a total 
of more than 5200 spectra of 169 systems, with from 5 to 



55 observations per target. Of these 169 systems, ~40% do 
not have observations at or near quadrature phases (even in 
some cases with 55 spectra!). The remaining ~100 stars were 
studied in greater detail, using the procedures described in 
Section [jj] 



3 DATA ANALYSIS 

The spectra were normalized by using an interactively de- 
fined spline fit to line-free continuum regions. (In examining 
each spectrum, any remaining noise spikes caused by cosmic 
rays etc. could be removed.) The spectra were then velocity- 
corrected to the heliocentric reference frame, and 'windowed' 
to isolate specific regions which include strong absorption 
lines, such as 4300-4500A (H7, He 1 A4388, A4471). Or- 
bital phases were calculated from photometric ephemerides 1 
published by Wyrzykowski et al. (2004, based on 4 sea- 
sons of OGLE-II photometry), where available; otherwise 
we adopted results from Udalski et al. (1998; 1^ seasons). 

3.1 Orbital parameters and component spectra 

H 3 03 described in detail the implementation of spectral dis- 
entangling codes which permit precise and direct determi- 
nations of the orbital parameters of the binary components, 
without the need for intermediate steps (i.e., measurement 
of velocities for each component for every observation). We 
demonstrated that our codes, which implement the algo- 
rithm described by Simon & Sturm (1994), give results 
that agree well with orbital solutions based on the cross- 
correlation methods previously more frequently used for 
measuring radial velocities. Compared to these 'traditional' 
methods, spectral disentangling not only yields orbital solu- 
tions of better precision (or, indeed, orbital solutions from 
data that would not be amenable to cross-correlation analy- 
ses), but also gives disentangled average spectra for each bi- 
nary component (at the expense of not returning 'per obser- 
vation' velocity measurements). Some excellent review pa- 
pers on the successes of spectral disentangling procedures 
are given by Hadrava (2004), Hensberge & Pavlovski (2004), 
Holmgren (2004) and Ilijic (2004). 

3.1.1 Velocities 

For each system we searched for the best (least-squares) or- 
bital solution by conducting grid searches in orbital semi- 
amplitudes (Ki,K2), given the period and reference epoch. 
Circular orbits were assumed unless the light curve demon- 
strated eccentricity. 2 It was found that using data in the 
phase ranges 0.20-0.30 and 0.70-0.80 gave the most reliable 

1 In consequence, when we refer to the 'primary component' we 
mean the star which is at superior conjunction at photometric 
phase zero (i.e., eclipsed at primary eclipse), which is the compo- 
nent with greater surface brightness - not necessarily the more 
massive component. 

2 In most cases, the light curves show no evidence of departures 
from circular orbits, with secondary eclipses at phase 0.50, and 
both eclipses of equal durations. In a few cases (discussed in Sec- 
tion the orbits are, however, detectably eccentric, with sec- 
ondary eclipses displaced from phase 0.50 by up to 0.10. 
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and consistent results; as necessary, wider phase ranges were 
employed (0.15-0.35 and 0.65-0.85), but observations closer 
to conjunctions offer no useful leverage on orbital-velocity 
amplitudes. 

At our spectral resolution, many binaries are clearly 
double-lined at quadrature phases, and well-defined orbital 
solutions are then readily obtained. However, the disentan- 
gling procedure is somewhat sensitive to the combination 
of orbital phase coverage and S/N in the data - important 
limiting factors in our sample. Where the grid search yielded 
multiple x 2 minima, the consequent ambiguity could some- 
times be safely resolved by further investigation (identifying 
prefered solutions on physical grounds). In all cases, sequen- 
tially finer grid searches were used to refine the solution. 

The outcome of this survey was that, of the ~100 sys- 
tems examined in detail, 40 yielded reliable orbital param- 
eters (in addition to the 10 systems already discussed in 
H 3 03). The remainder did not provide sensible solutions, in 
that derived masses and distances were obviously discrepant 
(typically by factors of 2 or more). Spectroscopic solutions 
for the 40 'good' systems are given in Table Q for circular 
orbits, and Table [2] for eccentric orbits. As the disentan- 
gling procedure returns only the semi-amplitudes, the tab- 
ulated systemic velocities were determined separately, from 
a cross-correlation analysis of the disentangled mean spec- 
tra against a template spectrum of the 09 V star 10 Lac, 
using the VCROSS code (Hill 1982). The values reported for 
each binary are the flux-weighted means of the two com- 
ponents. For eccentric orbits, values of eccentricity, e, and 
longitude of periastron, w, were determined from the light- 
curve solutions (from the position of secondary eclipse and 
the durations of the two eclipses), and then held fixed in the 
spectroscopic solutions (because of selective phase coverage 
of the spectroscopy). 3 

The published velocity semi-amplitudes of the 10 sys- 
tems studied by H 3 03 were derived from spectra obtained 
only in the first, 2001, observing session. Since the num- 
bers of useable spectra of these 10 systems has subsequently 
increased (and the original data have been re-reduced), 
new orbital solutions were conducted. The mean differ- 
ence between the published and new semi-amplitudes is 
—2 ± 4 km s _1 (s.e.), and in no case did the original and 
new semi-amplitudes differ by more than the joint 1-er un- 
certainties. The mean difference in systemic velocities is 
-ll±4kms- 1 (s.e.,H 3 03 minus this paper). This difference 
reflects the fact that in H 3 03, the heliocentric velocity cor- 
rections (~ 6 km s _1 ) were inadvertently applied with the 
incorrect sign, such that the previously published 7 veloci- 
ties should be corrected by +12 km s _1 . It is pleasing that 
this known mistake should be recovered in the data anal- 
ysis, and that the expected uncertainties in the velocities 
(~ ±15 km s _1 ) should be consistent with the dispersion 
in the results. Other than the correction to the 7 velocities, 
there is clearly no need to replace the H 3 03 orbital solutions. 



3 The light-curve solutions require input (the mass ratio) from 
the spectroscopic orbit, and the spectroscopic solutions require in- 
put ('non-keplerian' corrections, as well as e and u) from the light- 
curve solutions. The velocity corrections rarely exceed 5 km s~ x 
and a single iteration between light-curve and spectroscopic solu- 
tions gave convergence (including any adjustments to e and ui). 



3.1.2 Spectral types 

The separated component spectra that result from the dis- 
entangling procedure have S/N values that are ~ \J n/2 
greater than for each of the n individual observed spectra 
used in the analysis. We constructed disentangled spectra 
for each system by using all observations with orbital phases 
0.15-0.35 and 0.65-0.85, and the full observed wavelength 
range. 

These disentangled component spectra, shown in FigsQ 
and|5| were classified by reference to the digital spectral at- 
las of Galactic 0,B stars by Walborn & Fitzpatrick (1990) 
and to Lennon's (1997) SMC B-supergiant study. This task 
was performed independently by all three authors in order to 
assess the uncertainties associated with the spectral types. 
Typically we found that our classifications agreed to within 
one temperature subtype, and the adopted modal values 
for the primary components are given in Table For the 
(normally weaker) spectra of the secondary components, we 
provide the range of spectral types assigned by the three 
authors. We found that the O-Bl stars were quite easy to 
classify, but that spectral classifications were more difficult 
for later B types due to the absence of He 11 lines and, 
in these low-metallicity stars, the lack of detectable Si 11 
and other metal lines (see the discussion by Evans et al. 
2004). However, our classifications of primary components 
all agree to within one subtype, and it is these classifica- 
tions that establish the mean effective temperatures of the 
primary components (and hence the system absolute mag- 
nitudes; Section^. In this sense the spectral types assigned 
to the secondary components are much less important, al- 
though they do serve as consistency checks for the effective 
temperatures of the secondary components. 



3.2 Light-curve analyses 

The light-curve synthesis code light2 (Hill 1979; Hill & 
Rucinski 1993) was used for all analyses of the OGLE-II 
light curves. 



3.2.1 Data sources 

Our initial modeling made use of J-band photometry pub- 
lished by Udalski et al. (1998), based on 1~ seasons 
of OGLE observations reduced using standard techniques 
(point-spread-function fits to stellar images). This 'OGLE- 
PSF' database afforded typically ~ 140 observations per 
light curve, with an average photometric error of 0.014 mag 
quoted by Udalski et al. (1998). The overall rms scatter of 
the observations about the model fits was 0.019 mag, rang- 
ing 0.011-0.028 mag for individual systems, but some of the 
OGLE-PSF light curves were not particularly well sampled, 
with some eclipses defined by as few as 10 observations. 

Wyrzykowski et al (2004) subsequently published re- 
sults from four seasons of OGLE /-band photometry, re- 
duced by using the 'difference image analysis' (DIA) tech- 
nique described by Zebrun et al. (2001). Not all our sam- 
ple stars are included in this OGLE-DIA database, but for 
those that are there is improved accuracy in the ephemerides 
(though most changes are quite small - ~ 1-2 parts in 
10° for Porb); the light curves are better sampled (typically 
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Figure 1. Disentangled component spectra for 20 binaries (results for the remainder of the sample are shown in Fig.|2j. The normalized 
spectra are plotted as relative flux, with the secondary offset by —0.5 continuum units for clarity. Wavelengths are in A; the strong 
absorption lines include He, H<5, H7, and He 1 AA4026, 4144, 4388, 4471, 4713, as well as He 11 AA4200, 4542, 4686 in some cases. 
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Figure 2. Disentangled component spectra for 20 binaries (results for the remainder of the sample are shown in Fig.Q. The normalized 
spectra are plotted as relative flux, with the secondary offset by —0.5 continuum units for clarity. Wavelengths are in A; the strong 
absorption lines include He, H<5, H7, and He 1 AA4026, 4f44, 4388, 4474, 47f3, as well as He 11 AA4200, 4542, 4686 in some cases. 
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Table 1. Circular-orbit solutions (including 'non-keplerian' corrections). Identifiers are from Udalski et al. (1998, PSF) and Wyrzykowski 
et al. (2004, DIA), and for each target there are n useful spectra, characterized by the stated S/N . 
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Table 2. Eccentric-orbit solutions; oj is the longitude of periastron measured from the ascending node for the primary component. 
Semi-amplitudes include 'non-keplerian' corrections. Identifiers are from Wyrzykowski et al. (2004, DIA) and Udalski et al. (1998, PSF). 
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~ 300 observations per target); and the re-reduced photome- 
try is more accurate (quoted average photometric error 0.011 
mag), at least for uncrowded sources. 

We therefore repeated the entire analysis for the 35 sys- 
tems in our sample that have DIA photometry, all of which 
are reported to be 'uncrowded' sources (Wyrzykowski et al. 



2004.) In general, we found only small changes in the system 
parameters, but improved accuracy; the overall rms scatter 
of the DIA observations about the model fits is 0.014 mag, 
ranging 0.005-0.025 mag for individual systems. (One ex- 
ception, 5 311566, is discussed in Section 13.3.171 ) The five 
systems for which we have only the PSF photometry still 
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have quite well-defined light curves, and satisfactory solu- 
tions. 

3.2.2 Methodology 

To carry out the light-curve solutions, mass ratios were fixed 
from the spectroscopic solutions, and standard values were 
used for limb-darkening coefficients (appropriate to the lo- 
cal temperatures at each grid point over the surfaces of 
both stars in each binary), for gravity-darkening exponents 
(for stars with radiative envelopes), and for albedos and 
electron-scattering fractions (from Hill & Hutchings 1971). 
Synchronous rotation was assumed in all cases. 

We solved each light curve for the two stellar radii, the 
orbital inclination, and the temperature of the secondary 
component (strictly the /-band flux ratio between the two 
stars, but characterized by black-body fluxes at the appro- 
priate temperatures and effective wavelengths). The tem- 
perature of the primary component was fixed at the value 
corresponding to its assigned spectral type, and the default 
starting assumption for each solution was that each system 
is composed of stars that are detached from their respec- 
tive Roche lobes. Semi-detached configurations, with the 
secondary star filling its Roche lobe, were adopted only when 
the detached solutions indicated a secondary-star volume of 
98-100% of its Roche lobe. 

The disentangled spectra offer an independent con- 
straint on the flux ratio: requiring the overall relative 
strengths of the components' line spectra to match expecta- 
tions based on spectral type gives an estimate of the contin- 
uum flux ratio (in the B band). Because this is a relatively 
blunt tool, and because the binaries studied in this paper are 
fainter than those discussed in H 3 03 (and consequently their 
spectra have lower S/N values), we did not impose these 
spectroscopic flux ratios on the light-curve solutions, but 
used them, where necessary, to distinguish between compet- 
ing, and sometimes statistically indistinguishable, solutions 
to the I-band light curves with detached or semi-detached 
configurations. 4 

Table El gives the adopted values for optimized parame- 
ters, and the observed and model light curves are shown in 
FigsEHSl together with the corresponding (O— C) plots. In 
90% of cases, the (O — C) curves show no systematic phase- 
dependent deviations, demonstrating that satisfactory fits 
have been obtained; the few exceptions are discussed indi- 
vidually below. We provide reduced chi-squared values, xi, 
in Table 01 as a quantitative measure of the quality of each 
fit to each set of data. 

3.3 Notes on individual binaries 

In this section we comment on various aspects of individual 
solutions. Targets are identified by entries in the Udalski 
et al. (1998) catalogue and, in parentheses, the OGLE DIA 
catalogue (Zebrun et al. 2001). 

Of the systems considered here, several have had light- 
curve solutions attempted by Clausen et al. (2003; 5 140701), 

4 For a given system, detached solutions almost invariably imply 
flux ratios that are larger than those of semi-detached solutions, 
since the differences in radii of the two stars are generally larger. 



by Wyithe & Wilson (2001, 2002; 4 103706, 4 163552, 
7 66175, 7 255621, 8 87175), or by Graczyk (2003; 4 103706, 
4 163552, 5 311566, 6 221543, 8 104222). Because their solu- 
tions were unconstrained by spectroscopy, we do not make 
detailed comparisons with, or comments on, their results, 
which must be regarded as substantially less secure than 
ours. 



3.3.1 1 99121 (DIA 003851.98- 733433.2) 

From 11 spectra in the wider phase ranges 0.15-0.35 and 
0.65-0.85, a single minimum was found in the disentangling 
solution. The light-curve solution converged immediately to 
a semi-detached configuration, with a photometric flux ra- 
tio of 1.2 confirmed by the relative strengths of the spectral 
lines. The model fits the eclipse curves well, but the ob- 
servational data just before first contact of primary eclipse 
are depressed relative to those just after fourth contact, and 
this asymmetry is reflected in the value of xt. Such depres- 
sions in light curves occurring just before the onset of pri- 
mary eclipse are observed quite frequently in semi-detached 
(Algol-type) binaries and are indicative of absorption by the 
gas stream from the inner Lagrangian point of the Roche- 
lobe-filling secondary star falling towards the primary star. 
Just before primary eclipse, the gas stream is seen in projec- 
tion across the face of the primary component, and causes 
the extra absorption. 

3.3.2 4 56804 (DIA 004633.14- 732217.0) 

Spectra in the wider phase ranges yielded multiple \ 2 min- 
ima in the disentangling solutions, none of which gave physi- 
cally sensible results, but the four spectra closest to quadra- 
tures gave a mass ratio (primary/secondary) q < 1.0, and 
the resultant light-curve solution provides an excellent fit 
to the data (xt = 0.73) and a B-band flux ratio of unity, 
consistent with the spectra. 

3.3.3 4 103706 (DIA 004725.50- 732716.7) 

A simple disentangling solution from spectra in the nar- 
rower phase ranges (0.20-0.30, 0.70-0.80) coupled with a 
detached solution to this well-defined light curve provides a 
flat (O — C) curve, rms scatter of 0.013 mag, and xl = 1-29. 
The light-curve solution predicts a B-band flux ratio of 1.8, 
consistent with the spectra. 

3.3.4 4 1104 09 (DIA 004700. 19- 731843.1) 

There were no spectroscopic observations very close to 
phases 0.25 or 0.75, but by using the complete spectra in 
the wider phase ranges a unique solution was found, at a 
mass ratio of 1.544. The light curve is asymmetric, with the 
first maximum being 0.03 mag brighter than the second, 
which limits the quality of the light-curve fit, as reflected 
in the reduced chi-squared of 2.80. Nevertheless, the shapes 
of both eclipses are well modeled, providing sensible stellar 
radii for this semi-detached system. The implied photomet- 
ric flux ratio of 0.6 agrees very well with the spectroscopic 
value. 
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Figure 3. /-band light curves from OGLE-II (data points), plotted over 1.5 cycles, and model fits from light2 (solid lines) for 10 SMC 
binaries; note that the vertical scale varies between panels. The (O— C) scatter about the fit to each light curve is shown in the lower 
sub panel. 
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Figure 4. /-band light curves from OGLE-II (data points), plotted over 1.5 cycles, and model fits from light2 (solid lines) for 10 SMC 
binaries; note that the vertical scale varies between panels. The (O— C) scatter about the fit to each light curve is shown in the lower 
sub panel. 
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Figure 5. /-band light curves from OGLE-II (data points), plotted over 1.5 cycles, and model fits from light2 (solid lines) for 10 SMC 
binaries; note that the vertical scale varies between panels. The (O— C) scatter about the fit to each light curve is shown in the lower 
sub panel. 
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Figure 6. /-band light curves from OGLE-II (data points), plotted over 1.5 cycles, and model fits from light2 (solid lines) for 10 SMC 
binaries; note that the vertical scale varies between panels. The (O— C) scatter about the fit to each light curve is shown in the lower 
sub panel. 
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Table 3. Light-curve solutions. System types (column 2) are detached, semi-detached, or (in one case) in contact. Primary temperatures 
are based on spectral types; secondary temperatures characterize the /-band flux ratio (Section 13.21 . Radii are given in units of the 
semi-major axis of the relative orbit; for semi-detached systems these are tied to the mass ratio. 
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3.3.5 4 163552 (DIA 004753.24- 731556.5) 

An obviously double-lined system giving a simple disentan- 
gling solution and a good detached solution to the well- 
defined light curve, yielding an rms scatter of 0.011 mag and 
xt = 0.82. The photometric flux ratio of 1.0 agrees with the 
spectra. 



3.3.6 5 26631 (DIA 004859.84- 731328.8) 

Disentangling solutions of the 10 spectra in the narrower 
phase ranges resulted in several x 2 minima, but only one 
(with a mass ratio of unity) gave physically sensible results 
when combined with the light-curve analysis. Statistically 
indistinguishable solutions to the light curve were obtained 
with a detached configuration (the primary being the larger 
star), and with a semi-detached configuration (with the sec- 



ondary being larger). In the detached case, the primary fills 
97% of its Roche lobe, the expected B-band flux ratio is 2.7, 
and the resultant astrophysical parameters are inconsistent 
with evolution models for single stars. In the semi-detached 
case, the secondary fills its Roche lobe and the B-band flux 
ratio of 1.4 is consistent with the spectra. The s-d solution 
was therefore adopted, giving a light-curve fit showing an 
rms scatter of 0.016 mag and = 0.92. 

3.3.7 5 60548 (DIA 004835.40- 725256.5) 

A simple disentangling solution was obtained from 22 spec- 
tra in the wider phase ranges. The light-curve solution con- 
verged to a semi-detached configuration with a photometric 
flux ratio of 1.7, consistent with the spectra. Although the 
rms scatter of the observations about the model fit is small, 
at 0.011 mag, the xi value of 1.92 reflects the small asymme- 
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try in the observed curve just before first contact of primary 
eclipse. 



3.3.8 5 95194 (DIA 004950.42- 731931.6) 

The two O-type stars are clearly seen in double-lined spec- 
tra, with a spectroscopic flux ratio of 0.8. The disentangling 
solution was obtained from 12 spectra in the wider phase 
ranges, mostly covering first quadrature, and demonstrate 
that the secondary component is the more massive star. The 
light-curve solution (xt = 0.39) was then determined with 
a secondary component larger than the primary but a little 
cooler, as indicated by primary eclipse being ~ 0.015 mag 
deeper than secondary eclipse. The resulting B-band flux 
ratio is consistent with that seen in the spectra. 



3.3.9 5 14 0701 (DIA 004943.08- 725109.0, HV 11284) 

A simple disentangling solution from 12 spectra in the wider 
phase ranges led to a detached configuration from the light- 
curve solution, and a photometric flux ratio of 2.3, consistent 
with the spectra. The rms scatter of 0.011 and xt °f 1-52 
indicate a good solution. 



3.3.10 5 180064 (DIA 005044.74-731739.9) 

The 15 spectra covering both quadratures in the narrower 
phase ranges provided a simple and well-defined disentan- 
gling solution, despite the substantial flux ratio (> 2). The 
light curve has a low amplitude (0.08 mag), but the model 
fit is excellent, with an rms scatter of 0.009 mag, and xt 
of 0.53. The predicted B-band flux ratio of 2.7 is consistent 
with the spectra. 



3.3.11 5 208049 (DIA 005045.00- 725844 4) 

All but one of the 16 spectra used in the first disentangling 
solution lie in the phase range 0.65-0.85, and several min- 
ima were found. The quoted result is from the spectra in 
the phase range 0.70-0.80, which provided the best light- 
curve solution, yielding a semi-detached configuration with 
a photometric flux ratio of 1.4 (consistent with the spectra). 
The quality of the fit is reasonable, with rms scatter of 0.016 
mag, and xt — 1-67. 



3.3.12 5 243188 (DIA 005118.78- 733015.8) 

The best disentangling solution was obtained with 18 spec- 
tra from the wider phase ranges of 0.15-0.35 and 0.65-0.85, 
and the large-amplitude light curve yielded alternative de- 
tached and semi-detached solutions with different flux ra- 
tios. The semi-detached solution requires a flux ratio of 0.85, 
which is consistent with the spectra (the secondary's spec- 
tral lines are slightly stronger than those of the primary). 
The quality of the fit to the light curve is excellent, with an 
rms scatter of 0.009 mag and xt = 1-01, and a small dip 
just before first contact of primary eclipse. 



3.3.13 5 255984 (DIA 005129.62- 732137.7) 

Spectra in the phase ranges 0.20-0.30 and 0.70-0.80 yielded 
a single sensible disentangling solution, and the low- 
amplitude (0.12-mag) light curve provides a good detached 
solution with a flat (O— C) curve, rms scatter of 0.011 mag, 
and xt = 0.81. The model flux ratio of 1.5 is consistent with 
the spectra. 

3.3.14 5 277080 (DIA 005111.68- 730520.3) 

A simple disentangling solution from 13 spectra in the nar- 
rower phase ranges, and a semi-detached configuration for 
the light-curve solution, provide the best match to the flux 
ratio of unity seen in the spectra. The (O — C) curve has some 
deviations around first contact of primary eclipse, reflected 
in the rms scatter of 0.016 mag and xt = 1-55. Nonetheless, 
the overall fit appears to be good. 

3.3.15 5 300549 (DIA 005123.57- 725224.1) 

There are indications of emission lines at H7, and for the dis- 
entangling solution we used only the He I lines, with phase 
ranges 0.20-0.30 and 0.70-0.80. The light-curve required a 
phase shift of —0.0050 before a good semi-detached config- 
uration was obtained, with a photometric flux ratio of 1.2, 
matching the line spectra. The rms scatter is 0.010 mag, 
with xl = 0.78. 

3.3.16 5 305884 (DIA 005120.17- 724942.9) 

The 17 spectra in the narrower phase ranges give one phys- 
ically sensible solution, and the light-curve analysis pro- 
vides a reasonable fit with rms scatter of 0.011 mag, and 
xt — 1.63. The predicted B-band flux ratio of 1.5 is consis- 
tent with the spectra. 

3.3.17 5 311566 (6 67289, DIA 005134.85- 724545.9) 

This system falls in the overlap of two OGLE fields, and so 
has two identifiers. In our initial analyses we solved the PSF 
light curves for '5 311566' and '6 67289' both separately and 
combined, obtaining closely similar results. 

The DIA light curve is very similar in appearance to 
the PSF results, clearly showing an eccentric orbit with deep 
eclipses, but unfortunately the number of observations defin- 
ing the shapes of both eclipses is very limited (Fig. 4g), and 
this factor constrains the quality of the overall solution. The 
light curve was solved as an obvious detached system at a 
mass ratio of unity for the usual set of parameters together 
with e,u. An orbital solution with these e,to values fixed 
provided the semi-amplitudes, and the light-curve solution 
was slightly revised. The paucity of observations through 
eclipses may contribute to the high rms of 0.036 mag and 
xt — 5.57, by far the worst in the entire sample. 

3.3.18 6 11141 (DIA 005203.96- 731849.3) 

A simple disentangling solution allowed a straightforward 
solution for the light curve, which has a rather low amplitude 
of 0.22 mag. The (O— C) curve is flat, with an rms scatter 
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of 0.009 mag and xt = 0.64. The photometric flux ratio of 
1.7 is consistent with the spectra. 



3.3.19 6 152981 (DIA 005241.88- 724622.4) 

Only the first-quadrature phases for this binary were ob- 
served, but an orbit solution was obtained straightforwardly. 
The light curve has deep eclipses, and a semi-detached con- 
figuration was required to obtain the B-band flux ratio of 
0.8 seen in the spectra. The rms scatter about the model fit 
is 0.014 mag, whilst xt = 1-75 reflects apparent slight asym- 
metries in the eclipse curves. A slight phase shift (~ 0.002) 
may improve the fit. 



3.3.20 6 180084 (DIA 005342.42- 732319.9) 

Only 8 spectra in the wider phase ranges were available for 
this system, but these gave a unique solution, indicating 
stars of similar masses. The well-defined light curve with 
deep eclipses is accurately modeled, giving an rms scatter of 
0.014 mag and xt — 0.79, and the photometric flux ratio of 
unity is confirmed by the spectra. 



3.3.21 6 221543 (DIA 005340.40- 725222.0) 

For the second system in the sample with an obviously eccen- 
tric orbit, the light curve was solved at an initial mass ratio 
of unity and including e,cu as free parameters. The spec- 
troscopic solution was obtained with these e, u) values fixed, 
from 22 spectra. The derived mass ratio of 1.025 clearly did 
not significantly change the subsequent re-solution of the 
light curve, which yields an rms scatter of 0.015 mag and 
xt = 0.91. The photometric flux ratio of 1.3 is consistent 
with the spectra. 



3. 3. 22 6 251 04 7 (DIA 005344 ■ 05- 723124 . 0) 

A simple disentangling solution was secured, and the sub- 
sequent light-curve solution proved to be straightforward, 
with a semi-detached configuration providing a photometric 
flux ratio of 1.1. The DIA light curve is free of the obvious 
contamination seen in the PSF results, although the derived 
parameters are similar. A phase shift of 0.0020 was applied 
to the DIA light curve to improve the fit to the light curve, 
which shows some anomalies in primary eclipse but a good 
fit through secondary eclipse. The rms scatter is 0.022 mag, 
and xt = 1.01. 



3.3.23 6 311225 (DIA 005402.00- 724221.6) 

A total of 20 spectra cover both quadratures well, but more 
than one solution was found. The adopted light-curve so- 
lution is for a semi-detached configuration with a B-band 
flux ratio of 1.4, consistent with the spectra. The rms scat- 
ter of 0.009 mag and xt — 0.89 indicate a good solution 
to this shallow light curve, which nevertheless has a small 
asymmetry in the descending branch of primary eclipse. 



3.3.24 6 319960 (DIA 005405.26- 723426.0) 

Only the second-quadrature phases were sampled, by 15 
spectra, but a solution was obtained straighforwardly. The 
spectra indicate a flux ratio of 0.4 and a semi-detached con- 
figuration. This result is confirmed in the light-curve solu- 
tion, which yields a reasonable fit, with rms scatter of 0.012 
mag and xt = 1-57. Again, there is a small dip just before 
first contact of primary eclipse. 

3.3.25 7 66175 (DIA 005438.22- 723206.2) 

A total of 20 spectra in the wider phase ranges, but mostly 
from second quadrature, were used to establish the orbital 
solution of this obviously double-lined system with a flux 
ratio less than unity. The light-curve solution gives a good 
match to the overall shape of the curve, but there is a signif- 
icant depression just before first contact of primary eclipse, 
which is reflected in xt — 4.79, although the rms scatter is 
tolerable at 0.016 mag. The light-curve solution for a semi- 
detached configuration correctly predicts the spectroscopic 
flux ratio. 

3.3.26 7 120044 (DIA 005531.57- 724307.8) 

An obvious double-lined system with well-sampled quadra- 
tures yielded a simple disentangling solution from 12 spec- 
tra. The light curve is also well modeled, by a model with a 
detached configuration, giving a flat (O — C) curve, rms scat- 
ter of 0.009 mag and xt = 0.81, and a flux ratio of unity, 
matching the spectroscopic value. 

3.3.27 7 142073 (DIA 005554 44- 722808.5) 

Various competing x 2 minima in the disentangling solution 
from 19 spectra in the wider phase ranges were reduced to 
one physically sensible solution with the use of the 6 spectra 
in the narrower phase ranges. The resulting mass ratio of 2.0 
provides a solution to the well-defined light curve with deep 
eclipses as a semi-detached configuration and a flux ratio of 
2.7, consistent with the spectra. The model fit is generally 
good, with rms scatter of 0.016 mag, but a xt as large as 
2.52 reflects the fact that the primary eclipse in the model 
is too deep by 0.05 mag. 

3.3.28 7 189660 (DIA 005637.30- 724143.3) 

A set of 10 spectra in the wider phase ranges yielded a sen- 
sible orbit solution. A phase shift of +0.0050 was required 
to improve the model fit to the light curve, which has deep 
eclipses. The overall rms scatter is 0.016 mag, and xt = 1-52 
reflects the fact that there is, again, a small dip in the ob- 
served curve just before first contact of primary eclipse in 
this semi-detached system. The photometric flux ratio of 1.3 
is consistent with the spectra. 

3. 3.29 7 1 93779 (DIA 005621 . 80- 723701 . 7) 

Both quadrature phases were well sampled, with 21 spectra, 
and the orbit solution was simple. The light curve is well 
defined, and the model fit to the data is very good, yielding 
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a flat (O — C) curve, an rms scatter of 0.012 mag and xt — 
0.94, and a flux ratio of 2.3, which is consistent with the 
spectra. 

3.3.30 7 255621 (8 30634) 

A second target that occurs in two OGLE fields, and the 
third system in this sample with an obviously eccentric or- 
bit, 7 255621 has a large orbital eccentricity (e = 0.198). 
As with the other eccentric-orbit systems in this sample, 
the number of observations defining the narrow eclipses is 
very limited, and this factor handicaps the accuracy of the 
solution. We solved the PSF light curves for 7 255621 and 
for 8 30634 separately, confirming agreement (particularly in 
the determinations of e, uj). We then combined the data into 
a single light curve with the 7 255621 ephemeris, since there 
were no obvious zero-point differences in the two datasets. 
The combined light curve was solved for e,w, and the re- 
sults adopted in the disentangling solution using 11 spectra 
in the wider phase ranges. The light curve was re-solved after 
the orbit solution to yield a good fit to the available data, 
with a flat O— C curve, an rms scatter of 0.021 mag, and 
xt — 1.64. The photometric flux ratio of 2.3 is consistent 
with the spectra. This system is not included in the DIA 
database. 

3.3.31 8 87175 (DIA 005830.98- 723913.8) 

Limited spectroscopic data gave 7 double-lined spectra in 
the narrower phase ranges to provide an orbital solution. 
The well-defined light curve, with deep eclipses, yielded a 
detached configuration, with a reasonably flat (O— C) curve, 
an rms scatter of 0.016 mag, and xt — 1-43, and a flux ratio 
of 1.1, in agreement with the spectra. 

3.3.32 8 1 04222 (DIA 005825.08- 721909.8) 

A total of 11 double- lined spectra in the wider phase ranges 
were available for the orbit solution. The light curve, with 
deep eclipses of nearly equal depth, has some inconsistent 
observations through both eclipses. The best solution, with 
a reasonably flat (O — C) curve, yields an rms scatter of 0.023 
mag, and xt = 2.03, reflecting the scatter in observations 
through the eclipses. 

3.3.33 8 209964 (DIA 010016.05- 721243.9) 

The obviously double-lined spectra sampled both quadra- 
tures with a total of 14 spectra in the wider phase ranges. 
The light-curve solution indicates a semi-detached configu- 
ration and gives a flat (O— C) curve, an rms scatter of 0.009 
mag, and xt = 1-38. The photometric flux ratio of 0.9 is 
confirmed by the spectra, and again there is a small dip in 
the observed light curve just before first contact of primary 
eclipse. 

3.3.34 9 1 0098 (DIA 010052.90-724748.6) 

The narrower phase ranges have 7 spectra, yielding an orbit 
solution with a flux ratio greater than unity. The light-curve 
solution requires a semi-detached configuration to match 



this flux ratio, giving a flat (O— C) curve, an rms scatter 
of 0.009 mag and xl = 0.74. 

3.3.35 9 47454 (DIA 010052.04- 720705.5) 

Six spectra observed at both quadratures provide the orbital 
solution. The light curve appears to be symmetric but the 
observations have substantial scatter, which is reflected in 
the rms scatter of 0.016 mag, and xt — 1-94. The photomet- 
ric flux ratio of 1.1, confirmed by the spectra, results from 
a semi-detached configuration; a detached solution was also 
found, but the implied photometric flux ratio of 4.6 is clearly 
ruled out by the spectra. 

3.3.36 9 64498 (DIA 010117.26- 724232.1) 

A total of 10 spectra in the wider phase ranges was needed 
to provide an orbital solution. The shallow light curve has 
an amplitude of only 0.17 mag, but is symmetric. The light- 
curve solution provides a reasonable fit with a mostly flat 
(O — C) curve, an rms scatter of 0.011 mag and xt = 0.80. 
The photometric flux ratio of 1.8 is consistent with the spec- 
tra. 



3.3.37 10 37156 

Fourteen weakly double-lined spectra in the wider phase 
ranges were used to find the orbital solution. The spectra 
show some emission in H7, but using only the He 1 lines 
gives essentially the same solution. The light curve is quite 
well defined, but has minimal coverage of the eclipse curves. 
Initial solutions indicated that secondary eclipse was slightly 
displaced from phase 0.50, and re-solutions including e,uj 
yielded a sensible small eccentricity of 0.013. The final so- 
lution provides a flat (O— C) curve, with an rms scatter of 
0.013 mag and xt = 1-06, and a flux ratio of 2.9 that is 
consistent with the spectra. 

3.3.38 1 94559 

Clearly double-lined spectra with a flux ratio of unity pro- 
vided an orbital solution from 8 observations in the narrower 
phase ranges. Although the PSF light curve has sparsely 
defined eclipses, the light-curve solution shows that a semi- 
detached configuration is necessary to provide a photometric 
flux ratio in agreement with the spectra. The rms scatter is 
large, at 0.029 mag, and xt = 3.96 confirms the limited 
quality of this star's light curve. 

3.3.39 10 108086 

Only 6 double-lined spectra, which indicate a flux ratio near 
unity, were available in the wider phase ranges to provide an 
orbital solution. The light curve requires a contact config- 
uration, and was solved by Rucinski's routines for contact 
binaries incorporated in light2, to yield a fill-out factor of 
F = 1.70 (as defined by Mochnacki 1984). A phase shift 
of +0.0040 was required to improve the model fit, result- 
ing in a flat (O— C) curve, rms scatter of 0.016 mag, and 
xt — 1-88. This deep-contact binary, comprising two O stars 
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with masses of 17 and 14M©, has an orbital period of only 
0.88 days - a remarkable system. 



3.3.40 10 110440 

The spectra show emission in the H lines, and the orbital so- 
lution (using the spectra in the wider phase ranges) is based 
only on the He 1 A4471 line. The light curve is symmetric 
and well defined, and the solution yields a flat (O — C) curve 
with an rms scatter of 0.011 mag and xt ~ 0-47 - one of the 
smallest in the sample. The photometric flux ratio of 1.5 is 
confirmed by the spectra. 



3.3.41 Other systems 

Of the remaining ~60 binaries similarly investigated, many 
failed to provide meaningful orbital solutions due to the spe- 
cific combinations of observed orbital phases and S/N values 
achieved in the spectra, as noted in Section Some spectra 
were significantly contaminated by emission, particularly in 
the Balmer lines, presumably from nebular sources. At the 
spectral resolution employed in these fibre-fed spectra, we 
were unable to remove these nebular lines unambiguously, 
and hence orbital solutions could not be obtained reliably. 
Nonetheless, there are evidently several intrinsically inter- 
esting binary systems amongst this set, notably 6 89617 and 
6 194856, both of which are composed of two O-type stars 
with orbital periods of, respectively, 0.87 day and 0.65 day. 
Clearly, such rare objects deserve to be studied more in- 
tensively in order to establish their properties and to test 
models for the formation of close binaries of high mass. 



4 ASTROPHYSICAL PARAMETERS AND 
EVOLUTIONARY CHARACTERISTICS 

Astrophysical parameters can be calculated directly from 
the orbital and the light-curve solutions. Results are pre- 
sented in Table E] for the components in detached systems, 
and in Tabled] for the semi-detached and contact systems. 
The 1-cr uncertainties on the derived masses, radii, surface 
gravities, temperatures and luminosities have been calcu- 
lated by the standard procedures of propagation of uncer- 
tainties (errors) from the data given in Table[T] TableEJ and 
TableOl (cf., e.g., Hilditch 2001). For all systems, an uncer- 
tainty of ±1000 K on the mean effective temperature of the 
primary componentwas adopted, based on the discussion of 
the spectral-type-effective-temperature calibration given in 
H 3 03. 



4.1 Detached systems 

In any projection of the Hertzsprung-Russell diagram, the 
coeval components of a pre-mass-exchange binary must lie 
on an isochrone. This expectation is borne out for our sample 
of SMC detached systems, as illustrated in Fig.0 where we 
compare our results for detached systems with the stellar- 
evolution sequences given by Girardi et al. (2000), choosing 
models for Z — 0.004 as appropriate to SMC metallicity, 
and the mass-surface-gravity plane for the comparison. The 



average mass ratio for the 21 binaries is 1.2, so that differ- 
ential evolution due to mass differences is not generally ex- 
pected to be important. There are a few exceptions, notably 
4 103706 (18+lOM©), 5 255984 (12+7M ), and 10 110440 
(II+6M0), for which the locations of the components are 
nonetheless consistent with isochrones to within the 1-cr un- 
certainties. 

A second, and more stringent, observational test of the 
models is to compare the empirical, dynamically-determined 
masses with those inferred from the evolutionary tracks in 
the HR diagram (Fig.0. The differences between the evolu- 
tionary mass and the observed mass for each of the 36 stars 
in detached systems are listed in Tabled in the sense evolu- 
tionary mass minus observed mass. The average difference is 
AM = —0.8 ±2.1 Mq (s.d.), and when combined with the 6 
stars in detached systems from H 3 03, AM = -0.9±2.2M Q 
(s.d.). Considering that the typical uncertainties on the ob- 
served masses are ~ ±1-2M Q , this agreement is excellent, 
and confirms that where discrepancies occur between 'evo- 
lutionary' and 'spectroscopic' masses (e.g., Herrero et al. 
1992), it is the spectroscopic results that are more question- 
able. 

A final, straightforward comparison can be made by ex- 
amining isochrones in the mass-luminosity plane (Fig. EJ; 
all stars should lie above the ZAMS. In practice, most of 
the observed components are distributed sensibly within 
their 1-cr uncertainties of the expectations from the mod- 
els. The obvious exceptions are the secondary components 
of 4 56804, 6 11141, and 10 37156, and the primary compo- 
nent of 4 103706. 

The overall conclusion from these comparisons is that 
agreement is remarkably good between observations and the 
Girardi et al. models, with only 4 exceptions in a sample of 
44 stars. 

4.2 Semi-detached and contact systems 

The distribution of the components of the 28 semi-detached 
systems in the mass-surface-gravity plane (Fig. HOjl could 
hardly be more different from those of the detached sys- 
tems. Note that the lines connecting the components of each 
binary are all at obtuse angles relative to the isochrones, 
in complete contrast to the detached systems. In the HR 
diagram (Figure [^5, the secondary components appear to 
be systematically older than the primaries on the basis of 
single-star models, and are all much larger than expected 
for unevolved stars of the same mass. This is a very clear 
signature that these stars are the mass-losing components 
seen in the post-rapid-mass-transfer phase of case-A evolu- 
tion (the longer- lived, slow-mass-transfer phase of case A). 
In our earlier discussion of orbit and light-curve solutions for 
individual binaries, it was noted that a substantial number of 
the semi-detached systems displayed clear evidence of small 
depressions (0.01 — 0.05 mag) in their light curves occurring 
just before the onset of primary eclipse. This evidence for 
absorption by a weak mass-transfer stream (elaborated in 
Section 13. 3. H occurs in ~ 40% of our semi-detached sam- 
ple, and may be taken as further evidence of these systems 
being in the slow phase of mass transfer. 

Whilst the locations of the observed secondary compo- 
nents of these semi-detached systems in Figure Hoi and Fig- 
urc II II clearly show that they are evolved mass-losing stars, 
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Figure 7. The primary and secondary components (filled and open symbols) of 21 SMC detached binary systems plotted with l-cr 
uncertainties in the mass-surface-gravity plane (18 from this paper, 3 from H 3 03). Lines join the primary and secondary of each binary. 
Also plotted are evolutionary tracks (solid lines) for single stars of masses 5, 7, 12, 15, and 2OM0 and isochrones (dotted lines) for ages 
of 0, 1, 2, 5, 10, and 20Myr, taken from the models of Girardi et at (2000) for Z = 0.004, appropriate for the SMC. 
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Figure 8. The primary and secondary components (filled and open symbols) of 21 SMC detached binary systems plotted with 1- 
a uncertainties in the log T ei j-log L/L@ plane (HR diagram; 18 from this paper, 3 from H 3 03). Also plotted are the zero-age main 
sequence and evolutionary tracks for single stars of masses 5, 7, 12, 15, 20, 30Mq, taken from the models of Girardi et al. (2000) for 
Z = 0.004, appropriate for the SMC. 
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Figure 9. The primary and secondary components (filled and open symbols) of 21 SMC detached binary systems plotted with 1- 
<r uncertainties in the mass-log L/Lq plane (18 from this paper, 3 from H 3 03). Also plotted are isochrones (solid lines) for ages of 
0,4, lOMyr from the models of Girardi et al. (2000) for Z = 0.004, appropriate for the SMC. 
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Figure 10. The primary and secondary components (filled and open symbols) of 28 SMC semi-detached systems and 1 SMC contact 
system plotted with 1-ct uncertainties in the mass-surface-gravity plane (22 from this paper, 7 from H 3 03). Lines join the primary and 
secondary of each system. Also plotted are evolutionary tracks (solid lines) for single stars of masses 5, 7, 12, 15, 20Mq and isochrones 
(dotted lines) for ages of 0, 1, 2, 5, 10, 20, 50Myr, taken from the models of Girardi et al. (2000) for Z = 0.004, appropriate for the SMC. 
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Figure 11. The primary and secondary components (filled and open symbols) of 28 SMC semi-detached systems and 1 SMC contact 
system plotted with 1-cr uncertainties in the log T ef j-log L/Lq plane (HR diagram) (22 from this paper, 7 from H 3 03). Also plotted are 
the zero-age main sequence and evolutionary tracks for single stars of masses 5, 7, 12, 15, 20, 30Mq, taken from the models of Girardi 
et al. (2000) for Z = 0.004, appropriate for the SMC. 
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Table 4. Astrophysical parameters for detached systems 
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-0.7 


7 
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12.5 ± 0.5 
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the observed primary components occupy the main-sequence 
band. A natural question is whether these stars, which are 
the mass gainers in case A evolution, have the observed val- 
ues of logT B ff and log L/Lq that are appropriate for single 
stars of the same (newly-acquired) mass. That is, does the 
evolutionary mass concur with the observed dynamical mass 
also for the mass-gaining stars? We have repeated the test 
carried out for the stars in detached systems on the primary 
components of the semi-detached systems, and the individ- 
ual differences are listed in Table |U The average difference 
from TableGOis AM = -1.0 + 4.4 M© (s.d.), which becomes 
AM — —0.1 ± 4.4 Mq (s.d.) when the 7 primaries in semi- 
detached systems from H 3 03 are included. 

This apparent agreement between evolutionary and ob- 
served masses would suggest that the mass gainers do settle 
into equilibrium states during the slow phase of case A with 
observed properties appropriate for their newly-acquired 
masses. The examination by Braun & Langer (1995), which 
also reviewed earlier work on accretion on to massive main- 
sequence stars, shows that the mass gainers may well adopt 
such characteristics, provided that their core regions become 
again well mixed with hydrogen and therefore have 'reju- 



venated' cores appropriate for their new total mass. How- 
ever, mass gainers may not always take that route, but 
rather develop chemical structures unlike those of single 
stars, where the core mass remains unchanged from before 
the mass-transfer process. Such stars would appear to be 
under-luminous for their newly-acquired mass. These obser- 
vational data would suggest that the mass gainers do have 
rejuvenated cores. However, the dispersion about AM is a 
factor ~ 2 greater than that for the detached systems, which 
suggests that the true answer could be more complex. 

There are two extreme outliers in the distribution, the 
primary components of 5 300549 (with 5M = -10.4M©) 
and 9 175323 (with 5M = +9.6 M ), and if they are re- 
moved from the sample AM = —0.1 ± 3.6 Mq (s.d.). [Sys- 
tematic errors in observed masses would lead to systematic 
errors in derived distances and these two stars are respec- 
tively at distance moduli of 18.52 and 18.84, neither of which 
are unusual for the observed distribution.] Perhaps the ex- 
tra dispersion in A M is due to some stars not having re- 
juvenated cores and therefore lower luminosities, but the 
uncertainties of the data set are not sufficiently small to re- 
solve it. Perhaps one or two of the larger positive S M values 
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Table 5. Astrophysical parameters for semi-detached and contact systems 
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could be due to binary systems being observed at the end 
of the rapid phase of mass transfer where evolution-model 
luminosities are larger than equilibrium values (Wellstein, 
Langer & Braun 2001). 

These results on SMC semi-detached binaries are sim- 
ilar to those for semi-detached systems in the Milky Way 
galaxy (MWG) compiled by Hilditch (2004). For that sam- 
ple of only 10 systems, the average difference is AM = 
+2 ± 5 M Q , although some of those discrepancies could be 
due to the adoption of different T e ff scales for O-type stars. 
This type of investigation should be pursued further, with 
more data of higher quality and on a larger sample, because 
it could help to resolve the still poorly known efficiency of 
semi-convective mixing in mass-accreting stars (Wellstein et 



al. 2001), and whether stellar cores do become rejuvenated 
in mass-gaining events. 

Regrettably, there are no published stellar-evolution 
models for mass-exchanging binary systems in cases A 
and AB at SMC metallicity (Z ~ 0.004). However, the 
work of Wellstein & Langer (1999), and Wellstein, Langer 
& Braun (2001) for solar-metallicity binaries (Z = 0.020) 
provides many detailed examples of what is expected for the 
evolution of close binaries with initial orbital periods of a few 
days and initial masses in our observed range. Examination 
of their tracks for conservative case-A evolution, followed 
by case AB at later stages, shows sensible concurrence with 
well-established semi-detached 0,B binaries in the MWG, 
with the observed secondaries identified as the mass-losing 
original primaries, and the observed primaries as the mass- 
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gaining original secondaries (cf. Hilditch 2004). The models 
indicate that initial orbital periods of order 2-5 days remain 
around those values during the early mass-reversal stages 
and into the slow phase of case A, which is expected to last 
of order 3-5 Myr before the case AB phase starts (when the 
orbital periods become much longer). 

Our SMC semi-detached binaries fit this overall pic- 
ture well. Our sample is defined by an upper orbital-period 
limit of P < 5 days, and a faint absolute-magnitude limit of 
Mv — —3, corresponding to a spectral type on the main 
sequence of about B2. Of our 50 binary systems, 21 are 
in detached states and 28 are in post-rapid-mass-transfer 
semi-detached states with properties that seem to be con- 
sistent with conservative case-A evolutionary models. The 
durations of the detached phases in the Wellstein & Langer 
and Wellstein, Langer & Braun models are ~3-12 Myr (de- 
pendent on the initial orbital period and masses). The du- 
rations of the slow phase of case A are ~l-5 Myr (again 
dependent on the particular initial conditions). Note that 
the semi-detached phases of evolution for some systems are 
as long as, or longer than, the detached phases of others. 

On the basis of these models, do we therefore expect 
an observed non-coeval sample to display an approximately 
3:4 division between detached and semi-detached systems? 
Are the particular characteristics of the stars consistent with 
conservative case-A evolution, or is there an overall demand 
for some mass loss from the binary systems to have occurred 
during the rapid mass-transfer phase that lasts only some 
few x 10 4 years? Binary-star evolutionary models at SMC 
metallicity are required to compare with observational re- 
sults on the substantial sample of eclipsing binaries we have 
presented. The only contact binary in the sample is in deep 
contact, with a Mochnacki (1984) fill-out factor F — 1.70, 
and will presumably merge into a single star of some 30M Q . 
Its origins are unclear, but it is not a unique object amongst 
SMC or MWG OB binaries. 



5 DISTANCES 

Our methodology for distance determinations is fully de- 
scribed in H 3 03. In brief, the photometric analysis gives the 
relative sizes of the components in each system, and the 
spectroscopic orbit puts these dimensions on an absolute 
scale. An effective temperature is assigned to the primary on 
the basis of its spectral type; its absolute magnitude the fol- 
lows from a light-curve synthesis code which utilizes surface 
brightnesses from line-blanketed model atmospheres (we use 
those given by Howarth & Lynas-Gray 1989). The light- 
curve solution gives the relative brightnesses of the compo- 
nents, and hence the absolute magnitude of the system. The 
true distance modulus is determined using /-band data, as 
the reddening is smaller, and the intrinsic fluxes less tem- 
perature sensitive, at longer wavelengths; extinction is de- 
termined from [B — I) colours. 

Because many of our stars are somewhat cooler than the 
(predominantly O-type) targets of Paper I, a single 'hot star' 
colour is no longer adequate; instead we adopted intrinsic 
colours from an approximate fit to model results, 

(B - I) = max[-0.6, 2.97 - 0.78 log(T off (l))]. 

(log(T e f{) > 4.3). In a handful of cases this leads to values for 



A{I) that are smaller than the expected foreground value of 
~ 0.09 (E(B — V) — 0.05); we nonetheless retain the com- 
puted values in order not to bias the mean distance modulus. 

Results are summarized in Table The mean true dis- 
tance modulus for this sample is 18.88 ±0.30 (s.d.), ranging 
18.29-19.48, and the mean extinction is A(I) = 0.179±0.070 
(E(B-V) = A(7)/1.8 = 0.099±0.040). The mean DM may 
be compared with the value of 18.89 ± 0.14 (range 18.71- 
19.13) obtained from 10 stars by H 3 03. 

In the present sample, there is no correlation of DM 
with A(I), m(I), or 7 velocity, but there is a modest cor- 
relation with M(I) (Pearson and Spearman correlation co- 
efficients indicating that the null hypothesis of no correla- 
tion can be rejected with ~99.5% confidence). Fig.ll2lshows 
that this correlation arises from a simple selection effect: 
the sample omits intrinsically faint systems with large de- 
rived DM, simply because such systems are excluded by our 
apparent-magnitude cutoff. This biases the mean DM to- 
wards a too-small value. Restricting the sample to M(I) 
brighter than —3.0 removes any correlation and gives a mean 
DM of 18.967 ± 0.311 from 29 stars, which we adopt as the 
best estimate from the present sample. 

Our results also show no statistically significant cor- 
relation of DM with celestial co-ordinates. However, Groe- 
newegen (2000) has recently analysed the intrinsic structure 
of the Magellanic Clouds, using Cepheid photometry and 
adopting a simple triaxial model. For the SMC, he finds a 
~2kpc variation in mean distance (~0.07 in DM) over the 
~2° spatial extent of our sample. Our results are consistent 
with this, a linear fit of distance along the projected 'axis' 
of the SMC suggesting that the NE extreme is ~0.1m closer 
than the SW. 

Groenewegen also discusses the line-of-sight depth of 
the SMC, finding an rms depth of 0.11m (3 kpc) for 
Cepheids. While the distance moduli from H 3 03 and the 
present sample are consistent with a single mean value, 
the present sample has a significantly larger dispersion 
(with 95% confidence). If the dispersion in the H 3 03 sam- 
ple (for which internal and external error estimates were in 
close agreement) were representative of the stochastic er- 
rors in the present sample, then the residual dispersion of 
~0.28m(~8kpc) in the present sample would reflect astro- 
physical effects. However, this is surely a substantial over- 
estimate of the true dispersion in our parent population, 
given the significantly larger observational uncertainties in 
the current dataset. We consider that our results are there- 
fore again consistent with the Cepheid studies in this re- 
spect. 

Taking the weighted mean of the two eclipsing-binary 
samples (H 03 and the present paper), our best estimate 
of the true distance modulus to the SMC is 18.912 ± 0.035 
(60.6 ± 1.0 kpc), where the quoted s.e. is, of course, an in- 
ternal error; the result is subject to systematic uncertainties 
at the ~0.1 mag level, as discussed in detail by H 3 03. 

Much work on the extragalactic distance scale is tied to 
the LMC distance modulus (e.g., the HST Hubble Constant 
Key Project; Preedman et al. 2001), because the structure of 
that galaxy has historically been regarded as better under- 
stood. We intend to address LMC binaries directly in future 
work, but the difference in distance moduli, SMC— LMC, is 
quite well documented as +0.5m (e.g., Udalski et al. 1999, 
Groenewegen 2000). Our SMC results therefore imply an 
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Figure 12. Distance moduli as a function of absolute magnitude (/ band). O- and B-type stars are identified, as are results from H 3 03. 
The diagonal line roughly corresponds to the magnitude cut-off in our sample, and illustrates a modest selection bias. 



Table 6. Distance determinations. M{I) and Iq are the absolute and observed /-band quadrature magnitudes, and DM is the true 
(reddening-free) distance modulus. The /-band extinction A(I) is discussed in Section l5l 
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LMC DM of 18.41 ±0.04 ±0.1, as compared with a value of 
18.50 adopted in the Key Project work. 



6 CONCLUSION 

Together with results from H 3 03, our sample of SMC bi- 
naries includes 50 systems. With only two or three obvious 
exceptions, the typical uncertainties are ±10% on masses, 
±4% on radii, and ±0.07 on log L - remarkably good values 
considering the limited spectral resolution and data quality 
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available for this study, and the use only of single-passband 
light curves. For comparison, in a recent review of the ob- 
served properties of high-mass eclipsing binaries, Hilditch 
(2004) adopted similar upper limits for the uncertainties on 
these quantities derived for Galactic 0,B binaries, and found 
only 22 systems that satisfied these criteria. The aims of this 
project have thus clearly been realised, in that we have pro- 
vided good determinations of fundamental parameters for 
the largest sample of 0,B type eclipsing binaries achieved 
for any galaxy (including the Milky Way), and have made a 
direct determination of the SMC distance modulus that has 
exceptionally high precision (and very good accuracy). We 
find good agreement between our observations and the evo- 
lutionary models, with no evidence for a 'mass discrepancy' 
in the models. 
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